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ABSTRACT 

We report the successful identification of the type of the supernova responsible for the supernova 
remnant SNR 0509-675 in the Large Magellanic Cloud (LMC) using Gemini spectra of surrounding 
light echoes. The ability to classify outbursts associated with centuries-old remnants provides a new 
window into several aspects of supernova research and is likely to be successful in providing new 
constraints on additional LMC supernovae as well as their historical counterparts in the Milky Way 
Galaxy (MWG). The combined spectrum of echo light from SNR 0509-675 shows broad emission and 
absorption lines consistent with a supernova (SN) spectrum. We create a spectral library consisting 
of 26 SNe la and 6 SN Ib/c that are time- integrated, dust-scattered by LMC dust, and reddened 
by the LMC and MWG. We fit these SN templates to the observed light echo spectrum using 
minimization as well as correlation techniques, and we find that overluminous 91T-like SNe la with 
Atoi5 < 0.9 match the observed spectrum best. 

Subject headings: ISM: individual(SNR 0509-67.5) — supernova:general — supernova remnants — 
Magellanic Clouds 



1. INTRODUCTION 

Over 100 years ago, a rapidly expanding nebula 
was photographed by Ritchey around Nova Persei 1901 
(jRitchcy 1901a^b, 1902), and it was interpreted as a light 
echo from the nova explosion (jKaptevnl Il902f ). Later 
modelling of the physics of the scattering and the geom- 
etry that lead to apparent superluminal e xpansion con- 
firmed this interpretation (jCouderd I1939D . Since then, 
light echoes (whereby we mean a simple scattering echo 
rather than fluorescence or dust re-radiati on) have been 
seen m the Galactic Nova Sagittarii 1936 (ISwopelll940D 
and the eruptive variable V838 Monocerotis ([Bond et al.l 



|2003[ ). Echoes have also been observed from extragalac- 
tic s upernovae, w i th SN 1987A being the most famous 
case (|Crottslll988l: iSuntzeff et al]|1988l: iNewman fc Rest! 
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20061) , but also iiicludin g SNe 19 91T (ISchmidt et all 
1991 'Spa rks et afl Il999f). 1993J_JSus:erman feCrottl 
2002: Liu et al.ll2003D. 19 95E (iQuinn'e t al. 200%1998bu 
ICappcUaro ct al. 2001') 2002hh dW clch 2007), and 
2003gd [Sugerman.2005i : iVan Dvk et al... 2006. 1. 

By simple scaling argumen ts based on the visib ility o f 
Nova Persei (jShklovskiil 19641: Ivan den Berghll965lll975D . 
light echoes from supernovae as old as a few hundred to 
a thousand years can be detected, especially if the illu- 
minated dust has regions of high density (> 10~^cm~'^). 
More sophisticat ed models of scattered light echoes have 
been published (|Chevalieil Il986t ISugermanI 120031: iPatatI 
l2005f ) but the tabulations do not predict late-time light 
echo surface brightness. 

The few targeted surveys f or light echoes fro m super- 
nova e (Ivan den Berghl Il966l: iBoffi et al.lll999D and no- 
vae ()van den Berghl 119771 : ISchaefeil 119881 ) have bee n un- 
successful. However, these surveys did not use digi- 
tal image subtraction techniques to remove the dense 
stellar and galactic bac kgrounds. E ven the bright 
echoes near SN 1987A (ISuntzeff et al.l IT988') at U « 
21.3 mag arcsec"^ are hard to detect relative to the 
dense stellar background of the Large Magellanic Cloud 
(LMC). 

During the five observing seasons allocated to Super- 
MACHO Project observations, the LMC was observed 
repeatedly using the Mosaic imager at the Cerro Tololo 
Interamerican Observatory (CTIO) Blanco 4m telescope. 
An automated image reduction pipeline performed high- 
precision difference-imaging from September 2001 to De- 
cember 2005. We discovered light echo systems associ- 
ated with three ancient SNe in the LMC. The echo mo- 
tions trace back to three of the youngest supernova rem- 
nants (SNRs) in the LMC: SNR 0519-69.0, SNR 0509- 
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67.5, and SNR 0509-68.7 (N103B). These three remnants 
have also been identified as T ype la events, base d on 
the X-ray spectral abundances ([Hughes et al.1119951 ). We 
have dated these echoes to events 400- 800 years ago usin g 
their position and apparent motion (Rest et al."2 005bl ). 
Such light echo systems provide the extraordinary op- 
portunity to study the spectrum of the light from SN ex- 
plosions that reached Earth hundreds of years ago, deter- 
mine their spectral types, and compare them to now well- 
developed remnant structures and elemental residues. 

We have obtained spectra of light echoes from each 
of the three light echo groups with the Gemini-South 
CMOS spectrograph. While the light echoes of SNR 
0519-69.0 and 0509-68.7 are in very crowded regions of 
the LMC bar, the light echo features associated with SNR 
0509-67.5 are in a much less confused area. They are also 
the brightest light echo features we have discovered to 
date. We have extracted the light echo spectrum asso- 
ciated with SNR 0509-67.5 applying standard reduction 
techniques. 

The stellar spectral LMC background cannot be com- 
pletely removed from the fainter light echo features of 
SNR 0519-69.0 and 0509-68.7. We have obtained multi- 
object spectra with using GMOS on Gemini South sep- 
arated in time by one year in order to subtract off the 
stationary stellar spectral background and retain the (ap- 
parently moving) supernova echo light. In this paper we 
discuss the analysis of the light echo spectrum of SNR 
0509-67.5, and we determine the SN spectral type. We 
defer the analysis of the light echo spectrum of SNR 0519- 
69.0 and 0509-68.7 to a future paper. 

2. OBSERVATIONS & REDUCTIONS 
2.1. Data Reduction of Imaging Observations 

The SuperMACHO Project microlensing survey mon- 
itored the central portion of the LMC with a cadence of 
every other night during the five fall observing seasons 
beginning with September 2001. The CTIO 4m Blanco 
telescope with its 8Kx8K MOSAIC imager and atmo- 
spheric dispersion corrector were used to cover a mo- 
saic of 68 pointings in an approximate rectangle 3.7 deg 
by 6.6 deg aligned with the LMC bar. The images are 
taken through our custom "VR" filter (Ac = 625 nm, 
SX = 220 nm, NOAO Code c6027) with exposure times 
of 60 s to 200 s, depending on the stellar densities. We 
used an automated pipeline to subtract PSF-matched 
template images from the most recently acquired image 
to search for variability (iRest et all l2005al ; iGarg et aj] 
[mFTtlMiknaitis et aUmKYK The resulting difference im- 
ages are remarkably clean of the (constant) stellar back- 
ground and are ideal for searching for variable objects. 
Our pipeline detects and catalogs the variable sources. 

While searching for microlensing events in the LMC, 
we detected groups of light echoes pointing back to three 
SNRs in the LMC, SNR 0519-69.0, SNR 0509-67.5, and 
SNR 0509-68.7 (N103B). The surface brightnesses of the 
light echoes ranged from 22.5 mag arcsec"^ to the detec- 
tion limit of the survey of about 24 mag arcsec"^. 

2.2. Spectroscopic Observations 

Images were obtained on UT 2005 September 7 in the 
r' band using the Gemini-S GMOS spectrograph cover- 
ing a 5.5 X 5.5 arcmin field centered on the brightest 



echoes associated with SNR 0509-675. These pre-images 
were used to design a focal plane mask which included 
slitlets on 9 echoes, 9 stars, and 28 blank sky regions. 
The slits were 1.0 arcsec wide. Spectroscopy was obtained 
using GMOS with the R400 grating, yielding a resolution 
of 0.8A and a spectral range of 4500-8500A. 

Spectroscopic observations were obtained on UT 
2005 Nov 7, 2005 Dec 6 and 2005 Dec 7. A to- 
tal of six hour-long observations were made. The 
data were taken using the no d-and-shuffle technique 
(|Glazebrook fc Bland-HawthornI i2001 ), with the tele- 
scope nodded between the on source position and a blank 
sky field located 4 degrees away (off the LMC) every 120 
seconds. In total, 3 hours were spent integrating on- 
source. CuAr spectral calibration images and GCAL fiat- 
fields were interspersed with the science observations. The 
nod-and-shuffle technique provides for the best possible 
sky subtraction despite strong fringing in the red for the 
CCDs. 

2.3. Data Reduction of Spectroscopic Observations 

The GMOS data were processed using IRAF" and the 
Gemini IRAF package. The GMOS data was written as 
multi-extension FITS files with three data extensions, 
one for each of the three CCDs in the instrument. We 
performed the initial processing by extension (that is, 
by CCD), waiting until after nod and shuffle subtrac- 
tion to mosaic the extensions into a single array. First, 
an overscan value was subtracted and unused portions 
of the array were trimmed. As the CCD was dithered 
after each nod and shuffle observation in order to reduce 
the effects of charge traps, a separate fiat field was ob- 
tained for each science observation. The fiat field images 
were fit with a low-order spline for normalization and 
then each science frame was flattened with the appropri- 
ate normalized flat. Both "A" and "B" nods in a given 
frame were flattened with the same flat. The Gemini 
IRAF package gnscombine was used to combine all the 
observations of a given light echo while also performing 
the subtraction of the nod and shuffle components. The 
individual extensions were then mosaicked into a single 
array. 

Since the slitlets were oriented at various position an- 
gles tangent to the bright portions of the echo, each indi- 
vidual slitlet had to be rectified using a geometric trans- 
formation derived from the CuAr calibration lamp spec- 
tra. The slitlets with the brightest echo features were ex- 
tracted by collapsing each slitlet along the spatial dimen- 
sion. Wavelength calibration provided by the CuAr cali- 
bration lamps were taken through the same mask. A low- 
order polynomial was fit to the calibration lamp spectra 
and the solution then applied to each slitlet. The three 
brightest individual slitlets could then be combined^*. 
A slitlet that had been purposely placed in an appar- 
ently blank portion of the field was then used to create 
a background spectrum to represent the diffuse spectral 
background of the LMC. This background was scaled and 

IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

Each slitlet was in a different physical location on the mask, 
so each spectrum has a slightly different central wavelength, neces- 
sitating wavelength calibration before combination. 
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Fig. 1. — Observed spectrum of a light echo assicated with 
SNR 0509-675 (cyan hne). We have removed the chip gaps (4838 A- 
4894A, 6235A-6315A), and various LMC background emission lines 
([O III], Ho, [N II], [S II]). The black line indicates the boxcar- 
smoothed spectrum. The Ijroad absorption and emission lines are 
consistent with SN spectra. 

subtracted from the one-dimensional combined spectrum 
of the echo. Despite these efforts, some residual spectral 
contamination remains, as evidenced by narrow emission 
lines in the final spectrum. 

We used the spectroph otometric standard LTT 4364 
(|Hamuv et al.lll992l Il994f ) to flux-calibrate the individ- 
ual spectra using our own routines in IDL. The stan- 
dard star was not observed on the same nights as the 
echo observations. The relative spectrophotometry is 
expected to be good at the 5% level based on our ex- 
tensive experience with the reduction of a large sam- 
ple of low-z Type la SNe for which w e have photom- 
etry and spectroscopv (|Mathesonl [2OO7I ). We also used 
LTT 4364 to remove telluric features from the spec - 
tra usin g tech n iques described by Wade fc Hornd ()1988( ) ; 
iBesselll (|1999[ ): iMatheson eFaLl ([20001). Figure [1] shows 
the reduced spectrum of the light echo associated with 
SNR 0509-675, which exhibits broad emission and ab- 
sorption features consistent with SN spectra. 

3. ANALYSIS 

Several authors have addressed the scattering of light 
off dust particles and its effect on the sur face brightness 
and spectrum of the resulting light echoes (Coudcrc 1939]: 
Chevalier 1986: Emmering fc Chevalie r 1989; Sugermajj 
20031: lPatatlF2005l:lPatat et al.ll2006l e^g.). The light echo 
spectrum can be described by the time-inte grated SN 
spectrum attenuated b y the scattering dust ijSugermanI 
120031: iPatat et al.1 120061) . First, we describe how we cre- 
ate our library of 26 time-integrated SN la spectra and 6 
time- integrated SN Ib/c spectra (see Section [3T|) . Since 
our goal is to compare SN spectra, we introduce and de- 
scribe methods to correlate and compare SN spectra in 
Section [321 and we test these methods by comparing the 
integrated SN la template spectra. We introduce a tech- 
nique that estimates the Atoi5 of a given SN la based 
on its correlation with other SNe la. In Section 13.31 we 
describe how the spectrum is attenuated by the scatter- 
ing dust. We then test in Section 13.41 our method to 
estimate the Ami5 of a SNe by correlating its light echo 
with dust-scattered, time-integrated SNe la template on 
the example of SN 1998bu. 

3.1. Time-integrated Spectra 



The observed light echo is derived from scattering of 
the supernova light off of dust sheets. These dust sheets 
have light time travel dimensions which are significant 
with respect to the duration of a supernova's pre-nebular 
phase. Thus, the light echo represents a time-integration 
of the supernova flux attenuated b y the scattering dust 
([Sugermanll200llPatat et al.ll2006f) . We create 26 time- 
integrated SN la spectra using the lightcu rve and spec- 
tral library of Matheson] ([20071 ): I Jha et al.l ([2006( 1. and 6 
time-integrated SN Ib/c using the data references shown 
in Column (10) in Table [2l We cannot use a simple 
integration algorithm, like the trapezoidal rule, since 
there are often significant gaps of sometimes up to 15 
days in coverage due to bright time, weather conditions, 
and schedule constraints. Simple interpolation over non- 
homogenous coverage would not correctly account for 
the non-linear shape of the SNe lightcurves. Thus we 
calculate for each day a spectrum as the weighted aver- 
age of its two closest in time input spectra and scale it 
so that the spectrum convolved with the V filter agrees 
with the magnitude from the lightcurve fit. Then we use 
the trapez rule to integrate. All these steps are simple 
weighted summations of the input spectra, and can be 
combined to a weighted sum of all the input spectrum. 
Since we have normalized the input spectrum so that the 
spectra convolved with the V filter have the same refer- 
ence magnitude V = 15, these weights have the added 
benefit that they indicate the contribution of each input 
spectrum to the final integrated spectrum. This gives 
us the means to test if one or two input spectra domi- 
nate the integrated spectrum due to imperfect coverage: 
If the maximum weight is large, it is an indication of 
such problems. We use this maximum weight as a tool 
to grade the integrated spectra. In detail, we perform 
the following steps on each SN: 

• The spectral templates Si at epochs ti are normal- 
ized so that the spectra convolved with the V filter 
have the same reference magnitude = 15. 

• We fit T/-band l ightcurve templates of SN la 
(Prie to et an i2006^ and SN Ib/c (P. Nugent online 
library -'^^) to the observed lightcurves. These fitted 
lightcurves range from -15 days to -t-85 days with 
respect to the B-band maximum. 

• For a given time tk, we find the spectra that are 
closest to this date in both time directions, t~^ and 

■ We estimate the spectrum as the time- weighted 
average of these two spectra with Cii{tk) = {t~2 — 

tk)/{it2-i7i) andcM = itk-i^) / For 
times before the first or after the last spectrum, we 
just use the first and last spectrum, respectively. 
Then we normalize the spectrum by n{tk) so that 
the spectrum convolved with the V filter agrees 
with the magnitude from the lightcurve fit. 

S[tk) = n{tk) * {cii{tk)Sii + ca{tk)Si2) (1) 

• We integrate the spectrum from -15 days to -f85 
days with respect to the B band maximum. Calcu- 
lating the integrated spectrum is thus just a linear 

http: //supernova. Ibl .gov/~nugent/nugent_templates .html 
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combination of the input template spectra Si 

F(A)=5]5(t,) (2) 



where the Wj are functions of n(ife) and Ci{t}~) 



(3) 



As outHned before, the dominant comphcation creat- 
ing the time-integrated spectra is that for a given SN, 
the epochs for which spectra are available are non- 
homogenous, and one or two spectra can end up dom- 
inating the integrated spectrum. If the maximum weight 
Wmax = maximum(wj) is large, it is an indication of 
such problems. To reflect the relative quality of this ef- 
fect, we grade our integrated spectra by requiring that 
for grade A, B, and C the maximum weight fulfills 
Wmax < 25%, Wmax < 35%, and Wmax < 45%, respec- 
tively. We find 13, 7, and 8 SNe la of grade A, B, and 
C, respectively (see Table [J). For the SNe of type Ib/c, 
we find 3, 2, and 1 of grade A, B, and C, respectively 
(see Table 121). 

3.2. Methods to Compare and Correlate Spectra 

Our ultimate goal is to find the time-integrated tem- 
plate spectra that match best the observed light echo 
spectrum in order to determine the (sub)type of the 
SNe. One possibility is to fit the template spectra to 
the observed light echo spectrum by a simple normaliza- 
tion and calculate the x^- The intrinsic problem with a 
X^-minimization fit of spectra is that already small but 
low spatial frequency errors such as those due to errors 
in dereddening or background subtraction can warp the 
spectrum and lead to a bad measure of fit. An alternative 
to the x^-minimization approach is the cross-correlation 
technique. In this paper, w e use an implem entation of 
the correlation techniques of lTonrv &: Davii (119791. the 
Supe r Nova IDentification code (SNIP: iBlondin fc Tonrvl 
[2006L [200l . We compare and correlate the template 
SN la spectra: Do the templates correlate better if they 
have similar Amis? Can we determine the Am 15 of a 
SNe template just by correlating it to the other SNe la 
templates? In the following we discuss the advantages 
and disadvantages of these techniques. 

3.2.1. ')(^ -minimization 

We compare each SN la template spectrum with all 
others by performing a x^-minimization fit with a nor- 
malization factor as the only free parameter. We use 
the variance spectrum as 1 — ct errors for the observed 
spectrum. Figure [2] shows the versus their difference 
in ATO15: Arni5(i) — Ami5(j) for all SN la with grade 
A ot: B. There is a trend that for smaller differences 
in Atoi5 the is smaller, indicating that indeed SN la 
with similar Amis have similar spectra. However, there 
is still a big spread. This is most likely due to errors in 
dereddening or galaxy subtraction introducing system- 
atic errors. Figur^ shows the spectra of SN 1999dq and 
SN 1999aa. The spectrum of SN 1999aa is very simi- 
lar to the spectrum of SN 1999dq, with the exception 
that the SN 1999aa spectrum is warped toward the blue. 
This difference is most likely not real but an artifact of 
dereddening or background subtraction. 




Fig. 2. — This figure shows the versus their difi'erence in 
Amis: Ami5(i) — AmisQ) for all SN la in category A and B. 
There is a trend that for smaller differences in Amis the is 
smaller, indicating that indeed SN la with similar Amis are more 
similar. However, there is still a large spread which is most likely 
due to errors in dereddening or galaxy subtraction. 




5000 5500 6000 6500 7000 7500 
A in Angstrom 

Fig. 3. — This figure shows the spectra of SN 1999dq (sohd blue 
line) and SN 1999aa (dotted red line). Note that the spectra of 
SN 1999aa and SN 1999dq have almost identical spectral features 
but that of SN 1999aa is warped toward the blue. This difference 
is most likely not real but an artifact of dereddening or background 
subtraction. 



In order to avoid these problems, one has to do a very 
careful reduction and analysis to obtain the template SNe 
spectra. 

• Interpolate the observed B and V light curves (no 
reddening corrections). 

• Warp the spectra to match the observed B — V 
color at a given epoch. 

• Calculate the K-correction from the warped spec- 
trum. 

• Apply calculated K-corrections and deredden the 
B and V light curves to get the intrinsic B — V 
color. 

• Finally correct the spectrum by reddening and 
warping it to match the intrinsic color calculated 
from the light curve. 

Similar techniques are being developed for SN la 
lightcurvc fitting to get distances to SNe la. We are cur- 
rently implementing the se techniqu e s and we will discuss 
their implementation in lRest et al.l (|2008f ). 
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3.2.2. Cross-Correlation of Spectra with SNID 

One alternative to the ^^"iiiininiization approach is 
the cross-correlation technique. In this paper, we 
us e an im plementation of the correlation techniques 
of iTonrv"^ Davis (197 8j), the Supe rNova IDentification 
code (SNID; Blond in fc Tonrvl l2006. '2007V In SNID, the 
input and template spectra are binned on a common log- 
arithmic wavelength scale, such that a redshift (1+z) cor- 
responds to a uniform linear shift in log A. The spectra 
are then "flattened" through division by a pseudo con- 
tinuum, such that the correlation only relies on the rela- 
tive shape and strength of spectral features, and is there- 
fore insensitive to spectral color information (including 
reddening uncertainties and flux mis-calibrations). The 
pseudo-continuum is fitted as a 13-point cubic spline 
evenly spaced in log wavele ngth between 2500Aand 
lOOOOA. We refer the reader to lBlondin fc Tonrvl ()2007D 
where the SNID algorithm is described in full detail. 
Finally, the spectra are smoothed by applying a band- 
pass filter to remove low-frequency residuals (wavelength 
scale > 300A) left over from the pseitdo-continuum divi- 
sion and high-frequency noise (wavelength scale < 50A) 
components. The main motivation for applying such 
a filter lies in the physical nature of supernova spec- 
tra, which are dominated by broad spectral lines (~ 
100 - 150A) due to the large expansion velocities of the 
ejecta (~ 10000 km s~^). A more detailed explanation 
of the spectrum pre-processing arid cros s-correlation in 
SNID is given bv lBlondin fc Tonrvl ()2007f ). 

The input spectrum is correlated in turn with each 
template spectrum. The redshift, z, is usually a free 
parameter in SNID (indeed, this code was developed 
to determine the redshift of high-z SN la spectra; 
iMatheson et al.ll2005t [Miknaitis et al.ll2007D . but can be 
fixed when the redshift is known, as is the case here for 
the LMC (z « 0.001). The quality of a correlation is 
determined by the rlap quality parameter, which is the 
product of the Tonry & Davis correlation height-noise 
ratio (r) and the spectrum overlap parameter {lap), r 
is defined as the ratio of the height of the heighest peak 
in the correlation function to the root-mean-square of its 
antisymmetric component, while lap is a measure of the 
overlap in restframe wavel ength space between the input 
and template spectra (see iBlondin fc Tonrvl [2007i ) . For 
an input spectrum with the restframe wavelength range 
[Aq, Ai], lap is in the range < lap < ln(Ai/Ao). In what 
follows, a "good" correlation cor responds to rlap > 5 
with lap > 0.4 ([M atheson ct all l2005l : iMiknaitis et all 
I2007t IBlondin fc 'Ibnrvii200&i2007l ). Note that these lim- 
its are not derived from the dataset in this paper, but 
from comparing single epoch spectra from low and high- 
z SNe. 

We correlate each pair of SN la spectra using the 
SNID terchniques in a similar fashion what we have done 
with the x^-niinimization. Figure 3] shows the rlap val- 
ues versus difference in Ami5 of the SNID correlated la 
template spectra. There is a clear trend toward stronger 
correlation for spectra with similar Amis. Note that 
most of the pairs with rlap > 10 have a difference in 
Amis of less than 0.25. 

The SNIDtechnique "flattens" the spectra before cor- 
relating them (see Section [3.2.2^ . The flattened spectra 
produced by SNID have the additional advantage that the 




Fig. 4. — This figure shows rlap versus their difference in Amis: 
Amis(j) — Ami5(j) for all SN la in category A and B. There is 
a clear trend toward stronger correlation for spectra with similar 
Amis. 
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Fig. 5. — This figure shows Xj versus their difference in Amis: 
Amis(j) — Amis(i) foi" ^-H SN la in category A and B. The cor- 
relation between the goodness of fit indicated by and Amis is 
significantly better than for x^- 

is more robust against errors in dereddening or back- 
ground subtraction. We calculate x} for all SN la pairs 
of flattened template spectra and plot it versus the differ- 
ence in Amis in Figure [5] The correlation between the 
goodness of fit indicated by x} a-nd Amis is significantly 
better than for x^. 

We use the template spectra to test if can be used to 
determine the Amis : For a given SN la template spectra, 
we calculate the xj with respect to all other SN la tem- 
plates. We find the template SNe with the smallest X/ 
(noted xj min)^ a'l'^ estimate the Amis by computing the 
error-weighted mean Amis for the three templates with 
the smallest X/ < 2 x Xy rain ■ Figure [6] shows the Amis 
determined by the lightcurve shape versus the Amis es- 
timated using the X/ (denoted as Ami5(xj). For SNe 
la with small Amis ^ l-l; the agreement between true 
and estimated Ami5(x^) is excellent. However, for SNe 
la with Amis > 1-1, there is a significant spread. The 
reason is that the sample of SN la with Amis < 1.1 com- 
prises both normal objects and overlumin ous ones with 
spectra similar to SN 1991T or SN 1999aa (| Jefferv et all 
119921 : [jha et al.ir2006 l). These latter objects have spectra 
that show large deviations from normal SNe la, especially 
around maximum light (where the impact on the light- 
echo spectrum is the greatest). At intermediate Amis 
(1.2 < Amis ^5 1-6), however, SN la spectra are simi- 
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FlG. 6. — This figure shows Amis determined by the lightcurve 
shape versus the Amis estimated using the denoted as 



Ami5(x)). The SN la subtypes "Ia-91T 



"la-norm" 



and "la- 

91bg" are indicated with filled red circles, open blue squares, and 
filled black diamonds, respectively. The cyan pentagon inidicates 
SN 1998bu. For SNe la with small Amis < 1-1, the agreement 
between true and estimated Amis(x^) is excellent. However, for 
^Ne la with Amis > l-li there is a significant spread. 

ar, and our approach does not enable a clear determi- 
nation. The two points with larg e Amis > 1-8 are both 
subluminous, 1 991bK-hke SNe la (|Garnavich et al.ll200l 
iJha et ai1l2006f l. The spectra show significant deviations 
from normal and overluminous SNe la around maximum 
light, yet the Amis values are systematically underesti- 
mated. The disagreement is a simple consequence of the 
lack of 1991bg-like SN la in our set of spectral templates. 
With only two such templates with Amis > 1-6 (see Ta- 
ble [U the error-weighted mean Amis of the three best- 
matching templates will systematically bias the Amis 
determination to lower values. In Figure [7] we show 
the difference between the true Amis and the estimated 
Amis(x/) versus Ami5(x/)- For Amis(x/) < l-l? the 
standard deviation of the estimated Amis(x^) compared 
to the true Amis is 0.05 magnitudes. The calculated un- 
certainties of Amis(x/) are slightly underestimated by a 
factor of 1.3. Therefore we consider Amis(x/) as a very 
good estimate of the true Amis for Amis(x/) ^ 1-1 with 
uncertainties smaller than 0.1 magnitudes. If the fitted 
Amis > 1-1, then it can only be said that it is unlikely to 
be a 1991T-like SN la (1999cl is the only 91T-hke SN la 
in our sample with Amis > 1.1; see Table [T]). More im- 
portantly, Figure [7] also shows that we are able to accu- 
rately determine the Amis for SNe la with Amis < l-l- 
Since all 1991T-like SN la fulfill this condition, we are in 
principle able to not only determine the Amis but also 
the SN la subtype for these objects. We will see in the 
next section that the light-echo spectrum presented here 
is most probably a 1991T-like SN la with dml5 < 0.9. 

3.3. Single Scattering Approximation 
Several authors have addressed the single s c atter- 



ing_ 



approximatio n for light e choes fCouderc' '1939'; 
Chev alier 1986; Emmering fc Ch evalier 1989; Sugermaii, 
2003; Patat 2 0051 e.g.). Following the derivation of 
Sugerman (|20d3[ ). the surface brightness Bsc of scat- 
tered light with wavelength A scattering at an angle 6 off 
dust at position r and thickness Az is 

BsciX, 0, r, Az) = F{X)nH{r)Gir, Az)S{X, 9) (4) 

Where F{X) is the time-integrated spectra of the SN, 
n// (r) is the number density of hydrogen nuclei, G(r, Az) 



S 0.1 

i °° 




Fig. 7. — This figure shows the difference between the true Amis 
and the estimated Amis(x^) versus Ami5(x^). The SN la sub- 
types "la-91T", "la-norm", and "la-91bg" are indicated with filled 
red circles, open blue squares, and filled black diamonds, respec- 
tively. The cyan pentagon inidicates SN f998bu. For Amis(x^) ^ 

1.1, the standard deviation of the estimated Amis(Xj) compared 
to the true Amis is 0.05 magnitudes. 

is a geometrical factor depending on the geometry be- 
tween the observer, SN, and dust, and the integrated 
scattering function S{X, 6) which is described in more 
detail in Section fS. 3. II 

Since for this paper we are only interested in the rel- 
ative fluxes, we can drop all terms that are not depen- 
dent on the wavelength, add reddening by the LMC and 
MWG, and the modelled, observed spectrum Fgcat bas 
the form 

FscatW ^ C„o™Ci^ir«(A)Ci;f ^(A)F(A)5(A, 0) (5) 

where C^j.^^ ^ and Cj^^^'~^ are the reddening by the 
MWG and the LMC, respectively, as discussed in Sec- 
tion [3?3i3l 

3.3.1. Dust Properties 

In order to get the total integrated scattering function 
S{X,9), we add up the integrated scattering function for 
each individual dust type: 



s{x,e) = Y,Sx{x,e) 



(6) 



X 



Following I Weingartner fc Draind (|2001f ). X can be s, cn, 
and ci for silicon dust grains, carbonaceous dust grains 
with neutral PAH component, and carbonaceous dust 
grains with an ionized PAH component, respectively, 
where PAH stands for Polycyclic Aromatic Hydrocarbon. 
For a given dust type X, the integrated scattering func- 
tion Sx is 



Sx{\0) 



hc,x{\a)ag^x{0,X,a)fx{a)da (7) 



where Qsc,x is the grain scattering efficiency, Gg — tto^ 
is the grain cross-section, /x(a) is the grain size distri- 
bution discussed in Section 13.3.21 and <i>x(^, A,a) is the 
iHenvev fc GreensteinI ()1941[ ) phase function 



^x{x,e) 



l-g\{X,a) 



(H-5i(A,a)-25^(A,a)cos0)3/2 



(8) 



with g{X,a) the degree of forward scattering for a 
given grain. We integrate Sx{X,0) for the individual 
grain types by using the extended Simpson's method 
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(jPress et alj I1992D . We interpolate the values for 
Qsc,x and gxi^,a) using the tables provided by B. 
T. Drained" jJD raine fc Lee' ^98^; 'Laor & Dr^Hi EMI 
IWeingartnerfc Drainc 2001; Li & Drainc 2001^ 

3.3.2. Dust Grain Size Distribution 

We use the models defined by IWeingartner &: Draind 
(|2001h which consist of mixture of carbonaceous grains 
and amorphous silicate grains. Carbonaceous grains are 
PAH-like when small (a < 10~'^miu), and g raphite-like 
when large (a > 10~^^m)7 L"i fc Draineir2001l) . The dust 
grain size distribution /(a) is written as 



with ngj.{a) is the number density of grains with size 
< a and uh is the n umber density of H nuc l ei in both 
atoms and molecules. IWeingartner fc Draind (|2001[ ) de- 
rives the size distributions for different line of sights 
toward the LMC. We adopt the parameters of their 
"LMC avg" model and bg = 2 x IQ-^ (Table 3, 
IWeingartner fc Draind l2001h . We can then calculate 
the size distributions for carbonaceous dust fci{a) ~ 
Cionfia) and fcn( a) = (1 - Cion).fia) using equa- 
tions 2, 4, and 6 of IWeingartner fc Draind {^01). The 
PAH /graphiti c grains are assume d to be 50% neutral and 
50% ionized ()Li fc Draind l2001h . thus C,o„ = 0.5. For 
amorphous silicate dust fs we use equations 5 and 6. 

3.3.3. Extinction 
The extinction can be expressed as 

logio Ce^tiX) = -OA{A{X)/A{V))RvE{B - V) (10) 

For the MWG, we calculate the extinction C^;^^^{X) set- 
ting Ry — 3.1 and Emwg{B — V) — 0.07, and by calcu- 
lating {A{\)/A{V)) using Equations 1-3 of lCardelli et all 
(|T989il . 

The average internal extinction of the LMC is E{B — 
V) — 0.1 (Bessell 1991), but different populations give 
different results. IZaritskvi (1999) finds a mean E{B — 
V) = 0.06 from red clump giants and E{B -V) ^ 0.14 
from OB types. They attribute this dependence to an 
age-dependent scale height: OB stars having a smaller 
scale height lie predominantly in the dusty disk. We 
use Elmc{B — V)= 0.05, half the average internal ex- 
tinction value, since the most likely position of the SN 
is somewhere halfway through the LMC. Since the light 
echoes are not in the Superbubble of the LMC, we use 
the average Rv.lmc = 3.41 value of the "LMC Average 
Sample" in Table 2. iGordon et all ([2003'). Even though 
the extinction curves in the LMC and SMC have similar- 
ities to the MWG extinction curves, there are significant 
differences. Thus w e calculate {A(\)IA(V )) for the LMC 
using Equation 5 in lGordon et al.l (|2003f) and the values 
in the "Average" row of the section "LM C Average Sam- 
ple" in Table 3 in 'Gor don eFall (|2003D . 

3.4. Testing the method with SN 1998bu 

Several 100 days after the explosion, the light curve 
of the type la supernova SN 1998bu suddenly flattened. 

http:/ /www. astro. princeton.edu/~drainc/dust/dust. did. html 



At the same time, the spectrum changed from the typ- 
ical nebular emission to a blue continuum with broad 
absorption and emission f eatures reminiscent of t he SN 
spectrum at early phases (jCappellaro et al.ll2001h . This 
was explained by the emergence of a light echo from a 
foreground dust cloud. A similar case is SN 1991T, but 
its light echo spectra are of significantly less signal-to- 
noise. 

We use SN 1998bu as a test cases to see if we can 
determine the Amis of the SNe la from their light echo 
spectrum. We correlate the template spectra with the 
light echo spectrum and estimate with our method the 
Ami5(x^) of SN 1998bu. We can then compare how close 
this estimated Amis is to the true, lightcurve measured 
Amis = 1.02 ± 0.04. This is the ultimate test if the 
method works on a real world example. 

For the light echo of SN 1998bu we assume that the 
reflecting dust is z = 70pc in front of the supernova, 
the host galaxy extir i ction is Ay"'* = 0.86, and Ry = 
3.1 (jCappellaro et al.ll200lD . Using these values, we can 
calculate time-integrated, dust scattered, and reddened 
template spectra by applying equation [5] for 28 SNe la, 
which we denote template spectra in what follows. As 
described in the previous Section r3. 2. 2[ we calculate the 
rlap and X/ values using SNID for the observed light echo 
spectra with the spectra templates. 

Figure [8] shows the rlap values of the SN 1998bu fight 
echo and the template spectra for the different SN la 
subtypes In general, 91T-like SNe la (filled red circles) 
have a lower than normal Amis and are overluminous, 
slow decliners, whereas 91bg-like SNe la (filled black di- 
amonds) have a large Amis and are under luminous, fast 
decliners. Besides one outlier at Amis = 1-7, all other 
rlap > 5 values are for SN la with Amis < 1-2. This 
is in very good agreement with the Amis = 1-02 ± 0.04 
of SN 1998bu determined from its lightcurves (Note that 
the limit rlap > 5 is not derived from the dataset in this 
paper, but from comparing single epoch spectra from 
low and high-z SNe). Similarly, the best X/ values are 
all for spectra templates with Amis < 1-2 (see Fig- 
ure |9|). We apply our method to determine the Amis 
described in Section l3 . 2 . 21 using the X/ and estimate the 
Amis(x/) = 1-01 ± 0.03 for SN 1998bu. This is within 
the errors to the Amis determined with the lightcurves. 
Figure H and [7] show the Ami5(x/) of SN 1998bu (red 
circle). We find that the SNID correlation technique pro- 
vides more than sufficient discrimination between input 
templates for Amis < 1-1 

4. DISCUSSION 

We have obtained a spectrum of the light echo lo- 
cated at RA=05:13:03.77 and DEC=-67:29:04.91 at 
epoch J2000. The associated SNR 0509-675 is 
at RA=05:09:31.922 and DEC=-67:31:17.12 at epoch 
J2000. The angular distance between the light echo and 
the SNR is 0.340 degrees. Using the SNR age of 400 years 
determ ined using light echo apparent motion ([Rest et al.l 
I2005hl ^. we can determine the line-of-sight distance be- 
tween the dust sheet and the SNR. We find this distance 
to be 2; = 300 pc, which we use to calculate the scatter- 
ing angle 9. Then we can calculate the time-integrated, 
dust-scattered, and reddened template spectra by apply- 
ing equation [5] for 28 SNe la and 6 SNe Ib/c. We use 
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FlG. 8. — Correlation parameter rlap of the light echo of 
SN 1998bu determined with SNID versus Amis of the time- 
integrated, dust-scattered, reddened, and flattened SN la spectra. 
The SN la subtypes "Ia-91T", "la-norm", and "Ia-91bg" are in- 
dicated with filled red circles, open blue squares, and filled black 
diamonds, respectively. 



5000 5500 6000 6500 7000 7500 
A in Angstrom 

Fig. 10. — Time-integrated, dust-scattered, and reddened SN II 
spectra fitted to the observed light echo spectrum (cyan line) asso- 
ciated with SNR 0509-675. The chip gaps (4838A-4894A, 6235A- 
6315A), and various LMC background emission lines ([O III], Ha, 
[N II], [S II]) are removed from the observed spectrum. The fea- 
tures do not correlate, and we conclude that the observed light 
echo spectrum is not a SN II spectrum. 



• la-91T 
□ la-norm 

• la-91bg 




Fig. 9. — versus Amis for th^ fit of the time-integrated, 
dust-scattered, reddened, and flattened SN la spectra to the ob- 
served light echo spectrum associated with SN 1998bu. The SN la 
subtypes "Ia-91T", "la-norm", and "Ia-91bg" are indicated with 
filled red circles, open blue squares, and filled black diamonds, re- 
spectively. 

SNID to calculate the rlap and xj using both the SNe 
la and SNe Ib/c templates (see Column (9) and (10) in 
Table [Band Column (8) and (9) in Tabled). 

The spectrum of the light echo associated with SNR 
0509-675 shows broad emission and absorption features 
consistent with SN spectra. The question is, what kind 
of SN is it? Figure ^ shows the fit of two SN II tem- 
plate spectra (Type IIP and IIL) to the echo spectrum. 
Both spec tra are created by us ing a spectral library by 
Nugent^i (jGilliland et al.l 119991 ). The IIP spectra tem 



plate is based mostly on the models se en in [ Baron et al 



(2004), and its lightcurves are based on lCappellaro et al 



([1997). There is no need for any sophisticated correla- 
tion or fitting technique to conclude that the observed 
light echo spectrum is not that of a Type II supernova. 

Following the procedures described in Sec- 
tion 13.31 and 13.11 we have created the SN template 
spectra for 28 SNe la and 6 SNe Ib/c, and correlate 
these template spectra with the light echo spectrum 
using the SNID correlation technique described in Sec- 
tion 13.2.21 For the purpose of this paper, which is to 
identify the spectral (sub)type of the SN explosion that 
created SNR 0509-675, the SNID correlation technique 

http: //supernova. Ibl .gov/~nugent/nugent_templates .html 



described in Section 13.2.21 provides more than sufhcient 
discrimination between input templates. This technique 
"flattens" both spectra, and then correlates the main 
features of the spectra. The strength of the correlation is 
reflected in the parameter rlap, with values of rlap > 5.0 
indicating a strong correlation. Figure [TT] shows the rlap 
values versus Amis of the SNID correlated la template 
spectra. The SN la subtypes "Ia-91T", "la-norm", and 
"Ia-91bg" are indicated with filled red circles, open 
blue squares, and filled black diamonds, respectively. 
There is an improved match with template spectra 
with smaller Amis correlating more strongly with the 
observed light echo spectrum than the ones with large 
Atoi5. Only three templates, all with Amis < 0-9, 
have a strong correlation with the observed light echo 
spectrum. Figure [l2| shows a histogram of the rlap 
values for the different SN la as well as Ib/c subtypes. 
Note that the normal SNe Ic (SN 19941, SN 2004aw, 
SN 2005mf), denoted as "Ic-norm", show stronger 
correlations than the other types of Ib/c (the broad- line 
SN Ic SN 1997ef, the pecuhar SN lb SN 2005bf, and 
the normal SN lb SN 2005hg), denoted as "Ibc-other" 
in the lower panel of Figure [T^l (see also Table ^ . All 
SNe Ib/c have a significantly smaller rlap value than 
three of the 91T-like SN la SN, and no SN Ib/c has an 
rlap value bigger than 5, which is the cutoff value for a 
good correlation. 

The flattened spectra produced by SNID have the ad- 
ditional advantage that the is more robust against 
errors in dereddening or background subtraction. We 
calculate xj of the flattened template spectra with re- 
spect to the flattened observed light echo spectra for all 
templates (see Column (9) in Table [T] and Column (8) in 
Table [2|) . Figure [T3| shows xj versus Ami5 for the dif- 
ferent subtypes of SN la. The SN la subtypes "Ia-91T" , 
"la-norm", and "Ia-91bg" are indicated with filled red 
circles, open blue squares, and filled black diamonds, re- 
spectively. The correlation between xj a-nd Ami5 is, as 

expected, excellent: The 5 SNe with the smallest X/ also 
have the smallest Ami5. This confirms that rlap and x^ 
are equally suitable measures of fit when the low spatial 
frequency feautures are removed. Figure [T3| shows the 
histograms of xj for the different SN la subtypes (upper 
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FlG. 11. — Correlation parameter rlap determined with SNID 
versus Amis of the time-integrated, dust-scattered, reddened, 
and flattened SN la spectra. The SN la subtypes "Ia-91T", "la- 
norm", and "Ia-91bg" are indicated with flUed red circles, open 
blue squares, and filled black diamonds, respectively. An rlap 
value bigger than 5 is considered a strong correlation. There is 
a clear trend toward template spectra witli small Amis correlat- 
ing stronger with the observed light echo spectrum than the ones 
with large Amis. Only three templates, all with Amis < 0.9, 
have a strong correlation with the observed light echo spectrum 
(rlap > 5.0). 



1.0 1.2 1.4 1.6 1.8 

Ami-, 

Fig. 13. — versus Amis for the fit of the time-integrated, 
dust-scattered, reddened, and flattened SN la spectra to the ob- 
served light echo spectrum associated with SNR 0509-675. The 
SN la subtypes "Ia-91T", "la-norm", and "Ia-91bg" are indicated 
with filled red circles, open blue squares, and filled black diamonds, 
respectively. The 5 SNe with the smallest x^f a-r^ 91T-like SNe la 
and have the smallest Amis. Note that the goodness of fit cor- 
relates better with Amis than subtype. The 91T-like SNe with 
comparably large Amis have a larger Xf than the normal la with 
small Amis. 



Ia-91T 

la-norm 

la-91bg 



rlap (SNR 0.509-675) 

Fig. 12. — Histogram of rlap value determined with SNID for 
the different SN la subtypes ("Ia-91T", "la-norm", and "Ia-91bg", 
upper panel), normal SN Ic type ("Ic-norm", lower panel), and 
SN Ib/c type other than normal Ic ( "Ibc-other" , lower panel). An 
rlap value bigger than 5 indicates a good correlation between the 
two spectra. Note that the 91T-like SNe have the strongest corre- 
lation with the light echo spectrum. All SNe Ib/c have a signifi- 
cantly smaller rlap value than three of the 91T-like SNe la, and no 
SNe Ib/c has an rlap value bigger than 5, which is the cutoff value 
for a good correlation. 



panel) and SN Ib/c subtypes (lower panel). Similar to 
the rlap histograms, the normal SNe Ic have a decent 
xj, but all 91T-like SNe with Amis < 0.9 have a bet- 
ter Xf- We apply our method to determine the Amis 
of the light echo spectra using the and we find that 
Ami5(x/) = 0.87 ±0.05. 

Figure [15] shows the three 91T-like SN la template 
spectra with the best rlap values, overplotted on the ob- 
served light echo spectrum. The template spectra have 
the same features than the observed light echo spectrum, 
and that the agreement is very good. The only significant 
disagreement is that the template spectra have a slightly 
deeper absorption and stronger emission at 6IOOA and 
6500A, respectively. The two "normal" SNe la with the 
best rlap values are very similar and fit as well very good 
for wavelength smaller than 5800A, but the difference in 
strength of the two features at 6IOOA and 6500Ais more 



la-91T 

la-norm 

la-91bg 



Jtl 



I 



MJ 



□ Ic-norm 
■ Ibc-other 



I II 



X'j X 10 (SNR 0509-675) 

Fig. 14. — The upper panel shows the histogram of for the 
fit of the time-integrated, dust-scattered, reddened, and flattened 
SN la spectra fit to the flattened observed light echo spectrum 
associated with SNR 0509-675. Note that nearly all 91T-like SNe 
la have a better than any SNe of the other subtypes. The lower 

panel shows the same histogram of x^ for the SNe Ib/c. Note that 



the normal SNe Ic have a decent 
Amis < 0.9 have a better x^- 



X}, 



but all 91T-like SNe with 



pronounced (see Figure [TB)) . Note also that these two 
best-correlating normal SNe la are also the ones with the 
smallest Amis values. FigurefTTlshows the three SN Ib/c 
template spectra with the best rlap values. Even though 
there are similarities between the template spectra and 
the observed spectrum, the relative strength of the spec- 
tral features is not as consistent as for the SNe la. A 
few examples of SN template spectra with a poor corre- 
lation to the observed light echo spectrum are shown in 
Figure [m 

We conclude that the SN that created SNR 0509-675 is 
a 91T-like SN with Amis < 0.9. Knowing the SN type 
(and moreover its subtype and thus how energetic the 
explosion was) places stringent constraints on the explo- 
sion mechanism and hence on the interpretation of X-ray 
spectra of t he remnant. Analysi s of X-ray data of SNR 
0509-675 bv H ughes et al.l (|1995) classifies this SNR as a 
remnant of a SN la. Rec ent analysis of X-ray spectra by 
iBadenes fc Bravol (|2008i) also supports the classification 
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A in Angstrom 

Fig. 15. — The three time-integrated, dust-scattered, reddened, 
and flattened SN la spectra with the best rlap values (i.e. the best 
correlation with the observed light echo spectrum) , overplotted on 
the flattened observed light echo spectrum (solid cyan line). Note 
that all three template spectra are 91T-like SNe la with Amis < 
0.9. 




6000 6500 
A in Angstrom 



Fig. 16. — The two time-integrated, dust-scattered, reddened, 
and flattened spectra of normal SNe la with the best rlap values, 
overplotted on the flattened observed light echo spectrum (solid 
cyan line). Note that these two template spectra are also the two 
normal SNe la with the smallest Amis ■ 




6000 6500 
A in Angstrom 



Fig. 17. — The three time-integrated, dust-scattered, reddened, 
and flattened SN Ib/c spectra with the best rlap values, overplotted 
on the flattened observed light echo spectrum (solid cyan line). 
Note that all three template spectra are normal SNe Ic. 



as an over luminous, 91T-like SN la: models using hy- 
drodynamic calculations and nonequilibrium ionization 
simulations of highly energetic SN la reproduce the X- 
ray spectrum with its lin e flux ratios better than normal 
or subenergetic models (jBadenes fc Brav j |2008() . This 
is the first time that the (sub)type of an ancient SN has 
been determined by direct means by taking the spectrum 
of the original event. 




. light ectio 

snl998bp(!a-91bg) 
■ sn2005lig(lb-norm) 

snl997ef(lc-broad] 



5000 5500 



6000 6500 
A in Angstrom 



7000 



Fig. 18. — Three time-integrated, dust-scattered, reddened, 
and flattened SN template spectra of different types having a bad 
correlation with the observed light echo spectrum (cyan line). 

5. CONCLUSIONS 

We have obtained a spectrum of a light echo associ- 
ated with SNR 0509-675. By comparing and correlating 
time-integrated, dust-scattered, and reddened template 
spectra created from a spectral library of nearby super- 
novae of all types to the light echo spectrum, we find that 
overluminous, 91T-like SNe la match the observed spec- 
trum best. We correlate the template spectra with the 
observed spectra with SNID. The correlation parameter 
rlap is a measure of the strength of the correlation, and 
rlap > 5 indicates a strong correlation. Only SN la with 
Amis < 0.9 pass this cut. They correspond to intrin- 
sically overluminous SNe la wit h a spectrum res embling 
the prototypical SN 1991T ( Jef ferv et al1ll992f) . Simi- 
larly, all 91T-like SN la with Amis < 0.9 have a smaller 

than any other SNe when the "flattened" (see Sec- 
tion I3.2.2P SN templates are fitted to the observed spec- 
trum. Normal SNe Ic show some similarities to the ob- 
served light echo spectrum. However, the correlation is 
only weak (rlap < 5) and their is larger than the 

of the 91T-like SN la. Thus we can exclude them 
as the possible source event for SNR 0509-675. This is 
the first time that the (sub)type of a SN is conclusively 
and directly determined long after the event happened. 
Light echoes provide an excellent opportunity to connect 
the physics of the SN itself to its remnant. Much can be 
learned about the physics of SNe and their impact on the 
surrounding ISM from this direct comparison: Knowing 
the SN type (and moreover its subtype and thus how 
energetic the explosion was) places stringent constraints 
on the explosion mechanism and hence on the interpreta- 
tion of X-ray spectra of the remnant (jBadenes fc Bravol 
l2008f ). For the first time, models of SN explosions can 
now be tested for their mutual conistency with the SNe 
explosion itself and the observations of the SNR. We are 
currently working on expanding the sample of SNR with 
light echo spectra: In the LMC alone there are two more 
SNR with associated light echoes. Our investigation sug- 
gests also that subtyping of historical Milky Way super- 
novae, particularly the more recent SN 1572 (Tycho), 
SN 1604 (Kepler) and Gas A events, should be possi- 
ble provided that suitable light echo features are found 
and can be studied spectroscopically. Such a sample of 
SNRs with known explosion spectra will place stringent 
constrains on SNe explosion models and enhance our un- 
derstanding of these events that play such an important 
role in the production of heavy elements in our universe. 
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TABLE 1 



SN 
(1) 


Subtype 
(2) 


Amis 
(3) 




N 
(4) 


Pmin 

(5) 


Pmax 

(6) 


(7) 


Grade 
(8) 


Xf * w- 

(9) 


4 


rlap 
(10) 


snl991T 


Ia-91T 


0.94 ± 0, 


,05 


18 


-15.50 


71.50 


0.2041 


A 


2, 


,00 


1.8 


snl992A 


la- norm 


1.47 ± 0, 


.05 


11 


-7.50 


25.50 


0.2087 


A 


3, 


,74 


1.6 


snl997do 


la- norm 


0.99 ±0, 


.10 


12 


-12.05 


20.79 


0.2878 


B 


3, 


,98 


2.4 


snl998V 


la-norm 


1.06 ±0, 


,05 


8 


-1.99 


41.99 


0.3026 


B 


2, 


,55 


4.0 


snl998ab 


Ia-91T 


0.88 ±0, 


,17 


10 


-8.13 


46.74 


0.4090 


C 


1, 


,85 


2.3 


snl998aq 


la- norm 


1.05 ±0, 


,03 


14 


-1.24 


48.67 


0.2990 


B 


3, 


,42 


3.1 


snl998bp 


Ia-91bg 


1.83 ±0, 


,06 


11 


-5.03 


27.85 


0.2189 


A 


7, 


,37 


1.6 


snl998bu 


la-norm 


1.02 ±0, 


,04 


26 


-3.36 


43.66 


0.2159 


A 


3, 


,79 


3.4 


snl998dh 


la- norm 


1.23 ±0, 


,17 


10 


-10.07 


44.82 


0.4386 


C 


6, 


,16 


1.7 


snl998dm 


la-norm 


1.07 ± 0, 


,06 


10 


2.91 


45.86 


0.4007 


C 


3, 


,24 


4.2 


snl998ec 


la-norm 


1.08 ±0, 


,09 


6 


-3.98 


39.01 


0.2545 


B 


3, 


,88 


2.3 


snl998eg 


la-norm 


1.15 ±0, 


,09 


6 


-0.35 


23.68 


0.3581 


C 


3, 


,65 


2.7 


snl998es 


Ia-91T 


0.87 ±0, 


,08 


20 


-12.28 


43.72 


0.1753 


A 


1, 


,36 


5.4 


snl999aa 


Ia-91T 


0.85 ±0, 


,08 


22 


-11.16 


47.64 


0.2044 


A 


1, 


,32 


6.8 


snl999ac 


Ia-91T 


1.00 ± 0, 


,08 


16 


-5.97 


38.96 


0.1817 


A 


3, 


,16 


2.6 


snl999by 


Ia-91bg 


1.90 ±0, 


,09 


13 


-5.36 


40.67 


0.1568 


A 


6, 


,08 


2.0 


snl999cc 


la-norm 


1.46 ±0, 


,05 


7 


-4.14 


24.82 


0.2495 


A 


6, 


,17 


1.6 


snl999cl 


Ia-91T 


1.19 ±0, 


,19 


11 


-9.32 


36.67 


0.2792 


B 


2, 


,99 


1.5 


snl999dq 


Ia-91T 


0.88 ±0, 


,08 


21 


-13.02 


45.84 


0.1597 


A 


1, 


,41 


5.3 


snl999ej 


la-norm 


1.41 ±0, 


,05 


5 


-2.23 


10.76 


0.4068 


C 


4, 


,31 


2.4 


snl999gd 


la-norm 


1.16 ±0, 


,06 


5 


0.02 


33.89 


0.4036 


C 


5, 


,86 


1.9 


snl999gh 


Ia-91bg 


1.69 ±0, 


,05 


12 


2.02 


38.88 


0.4359 


C 


10, 


12 


2.5 


snl999gp 


Ia-91T 


0.87 ±0, 


,08 


8 


-6.29 


34.62 


0.1947 


A 


1, 


,64 


2.9 


sn2000cf 


la-norm 


1.27 ±0, 


,12 


6 


-0.16 


22.83 


0.3560 


C 


4, 


,23 


2.3 


sn2000cn 


Ia-91bg 


1.58 ±0, 


,12 


9 


-10.14 


26.77 


0.2764 


B 


6, 


,30 


2.6 


sn2000dk 


Ia-91bg 


1.57 ±0, 


,09 


6 


-5.11 


33.77 


0.2575 


B 


5, 


,07 


2.4 


sn2000fa 


la-norm 


1.00 ±0, 


,10 


13 


-12.03 


41.68 


0.1623 


A 


2, 


,42 


2.1 


sn2001V 


la- norm 


0.93 ±0, 


,05 


28 


-15.50 


48.50 


0.2020 


A 


2, 


,12 


4.4 



Note. — Overview of SN la template spectra jMatheson||2007l) . Column (1) shows the SN identifier. 
Column (2) shows the SN la subtype, where "la-91T" indicate the overluminous, slow decliners, "la- 
norm" arc the "normal" SN la, and "la-91bg" arc the undcrluminous, fast decliners. Column (3) shows 
the Amis of the SN la. The number of spectra N arc shown in Column (4), spanning a phase from 
Pmin to Pmin days (Columu (5) and (6)) with respect to the fitted V-band peak. Based on w^ax in 
Column (7) (The biggest weight assigned to a the spectra for a given SNc), a grade A, B, or C is 
assigned to each SNc (Column (8)), as described in Scction |3. 11 The Xf for the fit of the time-integrated, 
dust-scattered, reddened, and flattened SN la spectra to the observed light echo spectrum is shown in 
Column (9). Column (10) shows the rlap value determined with SNID indicating the correlation between 
the SN template and the observed light echo spectrum. An rlap value bigger than 5 is considered a good 
correlation. 



TABLE 2 



SN 

(1) 


Subtype 
(2) 


N 
(3) 


Pmin 

(4) 


pmax 

(5) 


^max 

(6) 


Grade Xf 
(7) 


* 10-* 
(8) 


rlap 
(9) 


Rcf. 
(10) 


snl994I 


Ic-norm 


18 


-6.50 


63.50 


0.2703 


B 


2.90 


3.8 


Filippenko ct al. (1995) 


snl997ef 


Ic-broad 


25 


-12.50 


82.50 


0.1485 


A 


4.08 


1.1 


jwamoto ct al. (2000) 


sn2004aw 


Ic-norm 


25 


-7.50 


46.50 


0.1610 


A 


2.17 


3.7 


Taubcnberecr et al. (2006) 


sn2005bf 


Ib-pec 


22 


-28.50 


33.50 


0.1474 


A 


6.68 


1.0 


Modiaz (2007) 


sn2005hg 


Ib-norm 


16 


-13.50 


26.50 


0.3285 


B 


7.45 


1.0 


Modiaz (2007) 


sn2005mf 


Ic-norm 


4 


3.50 


13.50 


0.4233 


C 


2.02 


4.2 


Modiaz (20d7) 


Note. — Same as Table [T] but with 
SN Ib/c, and Column (10) shows the re 


no Amis column, Column (2) 
ference of the SN data. 


shows the subtypes of 



